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ABSTRACT WhiA is a highly unusual transcriptional regulator related to a family of eukaryotic homing endonucleases. WhiA is 
required for sporulation in the filamentous bacterium Streptomyces, but WhiA homologues of unknown function are also found 
throughout the Gram-positive bacteria. To better understand the role of WhiA in Streptomyces development and its function as 
a transcription factor, we identified the WhiA regulon through a combination of chromatin immunoprecipitation-sequencing 
(ChlP-seq) and microarray transcriptional profiling, exploiting a new model organism for the genus, Streptomyces venezuelae, 
which sporulates in liquid culture. The regulon encompasses -240 transcription units, and WhiA appears to function almost 
equally as an activator and as a repressor. Bioinformatic analysis of the upstream regions of the complete regulon, combined 
with DNase I footprinting, identified a short but highly conserved asymmetric sequence, GACAC, associated with the majority of 
WhiA targets. Construction of a null mutant showed that whiA is required for the initiation of sporulation septation and chro- 
mosome segregation in S. venezuelae, and several genes encoding key proteins of the Streptomyces cell division machinery, such 
as ftsZ, ftsW, and ftsK, were found to be directly activated by WhiA during development. Several other genes encoding proteins 
with important roles in development were also identified as WhiA targets, including the sporulation-specific sigma factor o- whiG 
and the diguanylate cyclase CdgB. Cell division is tightly coordinated with the orderly arrest of apical growth in the sporogenic 
cell, and filP, encoding a key component of the polarisome that directs apical growth, is a direct target for WhiA-mediated re- 
pression during sporulation. 

IMPORTANCE Since the initial identification of the genetic loci required for Streptomyces development, all of the bid and whi de- 
velopmental master regulators have been cloned and characterized, and significant progress has been made toward understand- 
ing the cell biological processes that drive morphogenesis. A major challenge now is to connect the cell biological processes and 
the developmental master regulators by dissecting the regulatory networks that link the two. Studies of these regulatory net- 
works have been greatly facilitated by the recent introduction of Streptomyces venezuelae as a new model system for the genus, a 
species that sporulates in liquid culture. Taking advantage of S. venezuelae, we have characterized the regulon of genes directly 
under the control of one of these master regulators, WhiA. Our results implicate WhiA in the direct regulation of key steps in 
sporulation, including the cessation of aerial growth, the initiation of cell division, and chromosome segregation. 
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The filamentous, Gram-positive bacterium Streptomyces has a 
complex life cycle, culminating in the formation of specialized 
reproductive structures called aerial hyphae, which extend up- 
wards from the substrate mycelium into the air. Subsequently, the 
orderly arrest of aerial growth is coordinated with a remarkable 
synchronous septation event with concomitant segregation of 
chromosomes that leads to the formation of long chains of 50 to 
100 unigenomic spores (1-4). 

Genetic analysis of differentiation in the classical model organ- 
ism, Streptomyces coelicolor, defined two classes of developmental 
mutant, blocked at distinct stages of the life cycle, bid mutants are 
unable to form aerial hyphae and therefore lack the characteristic 
"fuzzy" appearance of the wild type (WT), exhibiting a shiny, 
"bald" phenotype (1, 3). In contrast, whi mutants are able to erect 



aerial hyphae but are unable to complete the life cycle to form 
mature chains of spores. Such mutants fail to synthesize the char- 
acteristic polyketide pigment associated with mature wild-type 
spores and therefore appear white on plates (1,3). 

The principal focus of this work is the WhiA protein. S. coeli- 
color whiA mutants fail to halt aerial growth or to initiate sporu- 
lation septation and the partitioning of chromosomes. Instead, 
they keep growing, forming extremely long, undivided aerial hy- 
phae containing uncondensed DNA that appears continuous on 
staining (5, 6). This phenotype is indistinguishable from that of 
whiB mutants, leading to suggestions that WhiA and WhiB func- 
tion together in the same regulatory pathway (5, 6). 

WhiA is of great interest, not only because of its critical biolog- 
ical role in controlling the initiation of sporulation but also be- 
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cause of the highly unusual nature of the protein itself, first pre- 
dicted bioinformatically (7) and dramatically confirmed when the 
structure of a WhiA homologue from Thermotoga maritima was 
determined (8). The N-terminal domain of WhiA, constituting 
two-thirds of the protein, is homologous to eukaryotic "LAGLI- 
DADG" homing endonucleases, but this domain is degenerate in 
that it has lost key residues required for metal binding and catal- 
ysis, and it also displays an extensively altered DNA-binding sur- 
face compared with true homing endonucleases (7, 8). Instead, 
DNA binding in WhiA is mediated through a smaller, C-terminal 
helix-turn-helix (HTH) domain that is completely absent from 
classical homing endonucleases (7-9). The two domains are con- 
nected by a flexible linker, and the C-terminal HTH domain binds 
DNA site-specifically by itself (9). Thus, the function of the much 
larger N-terminal degenerate homing endonuclease domain re- 
mains unknown. 

A further intriguing aspect of WhiA biology is that it is not 
confined to streptomycetes or even sporulating actinomycetes but 
is present in essentially all Gram-positive bacteria, including, 
therefore, many nonsporulating genera (for example, Staphylo- 
coccus). This observation, combined with its relationship with 
homing endonucleases, raises very interesting questions about the 
evolution of WhiA as a transcription factor (10) and about the 
reassortment of the target genes under its control in sporulating 
and nonsporulating bacteria. 

The regulons of the bid gene regulators BldA, BldD, BldH (also 
called AdpA), and, most recently, cr BldN have all been character- 
ized and reveal much of the network that governs the switch from 
vegetative to aerial growth (11-16). Thus, for example, we now 
know that cr BldN controls the expression of the chaplin and rodlin 
proteins, the principal components of the hydrophobic sheath 
that enables the aerial hyphae to escape surface tension and extend 
into the air (16). However, the whi gene regulatory network that 
controls the later stages of the Streptomyces life cycle remains 
largely uncharacterized. Consequently, our understanding of the 
regulation of processes such as the orderly cessation of aerial 
growth, the initiation of sporulation septation, and chromosome 
segregation remains incomplete. 

In-depth studies of the regulatory networks that control spo- 
rulation in Streptomyces have been greatly facilitated by the intro- 
duction of Streptomyces venezuelae as a new model system for the 
genus (16). The classical model system, S. coelicolor, poses certain 
logistical problems when seeking to characterize the regulons of 
the Bid and Whi master regulators. S. coelicolor sporulates only 
on solid medium, and the reproductive structures undergoing 
morphogenesis constitute just -5% of the total biomass, mak- 
ing the application of global techniques such as chromatin 
immunoprecipitation-sequencing (ChlP-seq) to development- 
specific transcription factors problematic. This study takes advan- 
tage of S. venezuelae, which, unlike S. coelicolor, sporulates syn- 
chronously in liquid culture to near-completion (17), greatly 
facilitating application of global techniques to the analysis of dif- 
ferentiation (16). We apply in vivo ChlP-seq and microarray tran- 
scriptional profiling to confirm the function of WhiA as a tran- 
scription factor and to characterize the regulon of genes directly 
under its control. These results implicate WhiA in the direct con- 
trol of key steps in sporulation, including the cessation of aerial 
growth, the initiation of cell division, and chromosome segrega- 
tion. 




FIG 1 Phenotypes of wild-type S. venezuelae (WT), the constructed AwhiA:: 
apr null mutant SV1 1 ( AwhiA::apr), SV1 1 carrying the empty vector pMS82 
{&whiA::apr affB <MT j::pMS82), and the complemented strain {&whiA::apr 
attB<p BTI ::whiA). Strains were grown on MYM solid medium and photo- 
graphed after 4 days. 



RESULTS AND DISCUSSION 

WhiA is required for sporulation in S. venezuelae. Before char- 
acterizing the WhiA regulon, a whiA null mutant of S. venezuelae 
was constructed to determine its phenotype and to compare it to 
the equivalent S. coelicolor whiA mutant (5). Using the Redirect 
PCR targeting method (18, 19), a whiA mutant was generated in 
which the coding region was replaced with an apramycin resis- 
tance (apr) cassette. As expected, the S. venezuelae AwhiA::apr 
strain was unable to sporulate on solid MYM (maltose-yeast 
extract-malt extract) medium or to synthesize the characteristic 
green spore pigment, instead appearing white (Fig. 1). In the wild 
type, the orderly arrest of aerial growth is followed by a septation 
event that divides the long, multigenomic hyphae into chains of 
unigenomic spores (Fig. 2A). To examine the effect of the whiA 
deletion in more detail, fluorescence microscopy of DAPI (4',6- 
diamidino-2-phenylindole)-stained cells and scanning electron 
microscopy were conducted. As seen previously for S. coelicolor 
whiA mutants, the S. venezuelae whiA mutant formed unusually 
long aerial hyphae that fail to initiate septation or chromosome 
segregation (Fig. 2B). It should be noted that wild-type S. venezu- 
elae forms straight aerial hyphae, in contrast to the coiled aerial 
hyphae made by wild-type S. coelicolor, and this is reflected in the 
phenotypes of the corresponding whiA mutants (Fig. 2B) (5, 6). 
All aspects of the S. venezuelae whiA null mutant were fully com- 
plemented in trans by a wild-type copy of whiA introduced into 
the <J>BT1 integration site (Fig. 1 and data not shown). Based on 
the phenotypes of the whiA mutants, we concluded that WhiA is 
likely to play similar roles in S. venezuelae and S. coelicolor. 

Defining the WhiA regulon. In order to determine which 
genes are directly regulated by WhiA in S. venezuelae, ChlP-seq 
experiments were conducted. We constructed a strain of S. ven- 
ezuelae that lacked whiA at its native locus but expressed an 
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FIG 2 Phenotypes of wild-type S. venezuelae (A) and the constructed AwhiA::apr null mutant SV11 (B), examined by scanning electron microscopy (i), 
transmission electron microscopy (ii), and fluorescence microscopy after staining membranes and DNA with 10 f^g/ml Nile red (iii) and 10 fig/ml DAPI (iv), 
respectively. Strains were grown on MYM solid medium for 2 days before fluorescence microscopy and 4 days before electron microscopy. 



N-terminal, triple FLAG-tagged version of WhiA from the 4>BT1 rulate in liquid culture. The best results were obtained when an 
integration site. This strain sporulated on solid MYM medium additional [Gly 4 Ser] 3 linker was engineered between the tag and 
(see Fig. SI in the supplemental material), although it displayed a WhiA. The resulting S. venezuelae strain sporulated to near- 
slight delay in spore pigment production and was slower to spo- completion in liquid culture, albeit slightly slower than the wild 
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type (Fig. SI and data not shown). Furthermore, 3XFLAG- 
[Gly 4 Ser] 3 -WhiA was readily detected using anti-FLAG antibody 
in Western blot assays and there was no visible cross-reaction with 
any proteins in wild-type S. venezuelae (Fig. S2), suggesting opti- 
mal conditions for clean ChlP-seq experiments. 

ChlP-seq was conducted with M2 anti-FLAG antibody at time 
points representing three distinct developmental stages: vegeta- 
tive growth, the onset of sporulation (fragmentation), and mid- 
sporulation. The wild-type strain (expressing nontagged WhiA 
from the native locus) was used as a negative control to eliminate 
any signals that might arise from cross-reaction of the antibody 
with other DNA-binding proteins. In addition, total (nonimmu- 
noprecipitated) input DNA from each time point was also sub- 
jected to sequencing. This very useful additional control enables 
nonuniform shearing of the chromosome to be taken into ac- 
count. More relaxed regions of DNA shear more readily than 
compact regions, and so the former can give rise to higher back- 
ground signals in ChlP-seq (20). Using P < 10~ 4 as the threshold 
for significance, a total of 284 peaks were detected in the FLAG- 
tagged WhiA strain over the three time points selected (see Fig. S3 
and Table SI in the supplemental material). Forty-eight of these 
were located more than 300 bp upstream of the nearest start codon 
(Table SI) and were not analyzed further. Using the same thresh- 
old for significance, only one peak was found in the wild-type 
(non-FLAG-tagged strain) control experiment. This peak (up- 
stream of sven5092) also appeared in the FLAG-tagged WhiA 
strain and was subsequently removed from the data set. 

Next, in order to determine how WhiA influences the expres- 
sion of its target genes, wild-type S. venezuelae and the congenic 
AwhiA mutant were subjected to time-resolved, genome-wide 
transcriptional profiling during vegetative growth and sporula- 
tion. Strains were grown under the same conditions used for the 
ChlP-seq experiments. RNA samples were prepared at 2-hour in- 
tervals from 8 to 20 h, by which time sporulation was nearing 
completion, and following cDNA synthesis and labeling, samples 
were hybridized to Affymetrix DNA microarrays. Three indepen- 
dent biological replicates were performed for each strain, and the 
transcriptional profiling data were examined for each of the genes 
potentially under WhiA control. Almost all the peaks identified by 
ChlP-seq were found in intergenic regions. Of these, 40% were 
found in intergenic regions between divergent genes. In the case of 
these divergent regulatory regions, the transcriptional profiling 
data set was used to assess which gene was a target for regulation by 
WhiA. 

Strikingly, 40% of WhiA target genes are downregulated in a 
whiA mutant and 35% are upregulated in a whiA mutant (see 
Table SI in the supplemental material). This strongly suggests that 
WhiA is bifunctional, working almost equally as an activator and 
as a repressor to control differentiation in S. venezuelae. 

Identification of a WhiA consensus binding site. To confirm 
and extend the ChlP-seq data and to help establish a consensus 
WhiA-binding sequence, DNase I footprinting analysis using pu- 
rified N-terminally histidine-tagged WhiA was carried out on ra- 
diolabeled probes derived from the upstream regions of six of the 
strongest WhiA targets, as judged by ChlP-seq. In each case, ad- 
dition of WhiA protein resulted in protection of a clearly defined 
region, consistent with high-affinity binding of monomeric WhiA 
(Fig. 3A). The sequences of these six confirmed WhiA target pro- 
moters were then examined using the MEME algorithm (21). This 
identified a short but perfectly conserved sequence (GACAC) in 



each of the submitted intergenic regions that correlated precisely 
with the position of the DNase I-protected region (Fig. 3B). This 
finding is in agreement with analysis conducted by Kaiser and 
Stoddard (9), in which WhiA binding was demonstrated at its own 
promoter. Following DNase I footprinting analysis, these authors 
tested the specificity of WhiA binding using a fluorescence com- 
petition assay and were able to determine those bases most impor- 
tant for the affinity of binding. Importantly, the GACAC motif 
independently identified here corresponds exactly with the se- 
quence identified by Kaiser and Stoddard (9). 

Following the identification of this short sequence, the remain- 
ing ChlP-seq targets were manually examined for the presence of 
the overrepresented WhiA-binding motif. Approximately 70% of 
the potential gene targets had a GACAC motif correlating closely 
with the position of the ChlP-seq peak. These target promoters 
were reexamined using the MEME algorithm, which failed to 
identify any additional overrepresented sequences associated with 
the GACAC motif (Fig. 3C). Within this subset, approximately 
50% of genes showed greater than a 2-fold change in expression 
when comparing wild-type S. venezuelae to the whiA mutant. 
However, WhiA ChlP-seq targets that lacked a GACAC motif 
were just as likely to exhibit clear whiA dependence. 

Because WhiA binds DNA as a monomer and it contains a 
large N-terminal homing endonuclease-like domain of unknown 
function, in addition to its small C-terminal HTH DNA-binding 
domain (8, 9), the orientation of the asymmetric GACAC motif 
relative to the target gene might well have regulatory significance. 
However, our analysis showed that there was no clear consistency 
in the directionality of the conserved motif, when present. We 
found that approximately twice the number of ChlP-seq- 
identified targets had upstream GACAC motifs compared to tar- 
gets with the reverse complement GTGTC sequence upstream. 
Out of those ChlP-seq targets that show clear whiA dependence in 
microarray analysis, -50% have upstream GACAC motifs, -25% 
have upstream GTGTC motifs, and -25% have no such sequence. 
Furthermore, this pattern does not change significantly when ac- 
tivated and repressed genes are analyzed separately. These data 
suggest that WhiA's ability to function as an activator or as a 
repressor is not tied to its orientation on the DNA. 

The identification of WhiA ChIP targets that depend on whiA 
for their expression, yet lack the GACAC motif, suggests that in 
some cases other sequences are sufficient to recruit WhiA. Three 
ChIP targets lacking the GACAC motif (sven3535, sven3229, and 
treZ, Table 1 ) were therefore tested for WhiA binding in vitro by 
DNase I footprinting. WhiA footprinted on only one of these three 
targets, treZ (sven5890). treZ is a clear WhiA target in ChlP-seq 
(Fig. 4A) that shows decreased expression in a whiA mutant, sug- 
gesting that WhiA activates transcription of this gene (see Table S 1 
in the supplemental material). DNase I footprinting of WhiA on 
the treZ promoter revealed two protected regions, both falling 
within the in vivo ChlP-seq peak (Fig. 4B). Although the determi- 
nants of this binding were not investigated, both protected regions 
contain the sequence GACGC, differing from the GACAC con- 
sensus sequence by 1 nucleotide (nt) (Fig. 4C). Thus, it is clear that 
WhiA can bind in vivo and in vitro to sites lacking the GACAC 
consensus sequence. For those WhiA ChlP-seq targets that did not 
bind WhiA in vitro, we speculate that additional transcription fac- 
tors or accessory proteins might be required. 

In Fig. 3, additional sites of weak protection were observed in 
the DNase I footprints on nrdR,ftsW, cslA, and thepaM-like gene 
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CAAACGGGCATACGCTCGATGACACGCCGCTTCGGCCCCTAGGATTCCCT 
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CCACTCGGACTACTCCAGTCGACACGCCGAGGGTTTGCCCCTTTTGTAGG 
GGTGAGCCTGATGAGGTCAGCTGTGCGGCTCCCAAACGGGGAAAACATCC 
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|1- 



f 



CMCO*in<ON00010 



FIG 3 DNase I footprinting analysis of WhiA binding to the promoter regions of nrdR,ftsW,filP, csIA, and svenl406 and the intergenic region between sven4724 
and sven4725. 5'-end-labeled probes were incubated with increasing concentrations of WhiA (indicated in ju,M above the lanes) and subjected to DNase I 
footprinting analysis as described in Materials and Methods. Footprints are flanked by Maxam and Gilbert sequence ladders (AG). Black brackets indicate the 
position of the WhiA-protected region. (B) Summary of DNase I footprinting results presented in panel A. Brackets indicate the protected sequences with the 
consensus GACAC motif shown in bold. The number indicates the distance (bp) between the consensus motif and the putative start codon of the downstream 
gene. (C) Overrepresented motif associated with -70% of WhiA target genes, identified using the MEME algorithm (21). The height of the letters in the sequence 
logo, in bits, is proportional to the frequency of the A, C, T, or G nucleotides at each position of the motif. Both the forward (i) and the reverse complement (ii) 
versions of the consensus are shown. 
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TABLE 1 Selected ChlP-identified WhiA targets 
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" Genes flanking the ChIP peak are listed. 

b Distance to the predicted start codon of the flanking genes. 

c (Possible) gene function based on annotation in StrepDB (http://strepdb.streptomyces.org.uk). 

d Maximum significance values (P < E— 04) for the 3 ChIP samples; V, vegetative growth; F, onset of sporulation/fragmentation; S, midsporulation. 



svenl406. It is not clear what these sites represent. Some contain 
GACAC motifs, and some contain the related sequence GACAG, 
but others do not show any obvious relationship to the consensus 
binding site. 

WhiA functions as an autoregulator. whiA appears to be the 
last gene in a four-gene operon, and previous work in S. coelicolor 
had suggested that whiA is transcribed throughout growth from a 
constitutive promoter located well upstream of the whiA gene 
itself (5). Consistent with this suggestion, our microarray tran- 
scriptional profiling data showed that S. venezuelae whiA is tran- 
scribed throughout vegetative growth and sporulation, and West- 
ern blot analysis of WhiA expressed from its native locus using 
polyclonal WhiA antiserum showed that the protein is similarly 
present throughout vegetative growth and sporulation (Fig. 5). In 
S. coelicolor, whiA expression is upregulated during development 
from a sporulation-specific promoter, whiAp2, lying just inside 
the upstream gene (5), and expression from whiAp2 is reduced to 
a low constitutive level in a whiA null mutant, suggesting positive 
autoregulation (22). Furthermore, the ability of WhiA to bind to 
the S. coelicolor whiAp2 promoter has been demonstrated in vitro 



(9). In line with these findings, our ChlP-seq experiments showed 
that whiA is a direct in vivo target of WhiA in S. venezuelae (Fig. 6A 
and Table 1). DNase I footprinting using purified N-terminally 
histidine-tagged S. venezuelae WhiA showed that WhiA binds to 
the S. venezuelae whiAp2 promoter region (Fig. 6B and C), in 
agreement with equivalent experiments carried out in S. coelicolor 
(9). Overall, it appears that WhiA binds to the region immediately 
upstream of the whiA gene in vivo to increase its expression prior 
to sporulation in Streptomyces. These and other experiments re- 
ported here suggest that WhiA activity is not a simple reflection of 
WhiA abundance, likely implying that WhiA activity is regulated 
posttranslationally (see "Conclusions"). 

WhiA targets involved in polar growth. Streptomyces hyphae 
grow by tip extension, and this polar growth is directed by a com- 
plex of coiled-coil cytoskeletal proteins called the polarisome (also 
called the tip organizing center or TIPOC). The three known po- 
larisome proteins are DivIVA, which is essential for viability of the 
organism, and two nonessential cytoskeletal proteins that never- 
theless play vital roles in tip growth, Scy and FilP (23-26) . All three 
members of the polarisome interact directly with each other, and 



6 mBio' mbio.asm.org 



September/October 2013 Volume 4 Issue 5 e00684-13 



Targets of WhiA in S. venezuelae before Sporulation 



A B treZ 



treZ AG 0 0.25 0.5 1.0 




forward 



FIG 4 treZ, a WhiA target lacking a perfect GACAC motif. (A) ChlP-seq data for the WhiA target gene treZ. Color coding of the ChIP samples is as follows: 
3XFLAG-[Gly 4 Ser] 3 -WhiA strain during vegetative growth (V, red), at the onset of sporulation/fragmentation (F, purple), and at midsporulation (S, blue) and 
the S. venezuelae negative-control strain (C, green). The plot spans approximately 3 kb of DNA sequence. Genes running left to right are shown in green, and 
genes running right to left are shown in red. The black arrow indicates the gene subject to WhiA regulation. (B) DNase I footprinting analysis of WhiA binding 
to the promoter region of treZ. 5' -end-labeled probes were incubated with increasing concentrations of WhiA (indicated in juM above the lanes) and subjected 
to DNase I footprinting analysis as described in Materials and Methods. The footprint is flanked by a Maxam and Gilbert sequence ladder (AG). Black and red 
brackets indicate the positions of the WhiA-protected regions. (C) Sequence summary of DNase I footprinting results presented in panel B. Black and red 
brackets indicate the protected sequences, with the near-consensus motif shown in bold. The number indicates the distance (bp) between the nearest consensus 
motif and the putative start codon of the treZ gene. 



Scy has been suggested to function as a molecular scaffold for the 
other components (26). The intermediate filament-like protein 
FilP is recruited by DivIVA and localizes just behind DivIVA at the 
hyphal tip (25). The reduced rigidity of the cell wall at hyphal tips 



in a filP mutant (revealed by atomic force microscopy) and its 
distorted morphology during growth suggest that FilP plays a me- 
chanical role in controlling cell shape (27). Consistent with this 
suggestion, FilP can form potentially load-bearing net-like struc- 



lOhr 12hr 14hr i6hr 18hr 

hA + A +A +A +A 

vegetative mycelium 



fragmentation 



WhiA 



sporulation 

FIG 5 Western blot analysis of WhiA abundance during submerged sporulation in wild- type S. venezuelae ( + ) and in the congenic whiA mutant control strain, 
SV11 (A). Strains were grown in MYM liquid sporulation medium, 10 jug total protein was loaded per lane, and the blot was developed using a polyclonal 
antiserum raised against S. venezuelae WhiA. Bars indicate vegetative growth, fragmentation, and sporulation in the wild type, as judged by microscopic 
examination. 
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forward 

FIG 6 WhiA functions as an autoregulator. (A) ChlP-seq data for WhiA binding to the whiAp2 promoter. Color coding of the ChIP samples is as follows: 
3XFLAG-[Gly 4 Ser] 3 -WhiA strain during vegetative growth (V, red), at the onset of sporulation/fragmentation (F, purple), and at midsporulation (S, blue) and 
the S. venezuelae negative-control strain (C, green). The plot spans approximately 3 kb of DNA sequence. The black arrow indicates the gene subject to WhiA 
regulation. (B) DNase I footprinting analysis of WhiA binding to the whiAp2 promoter. 5' -end-labeled probes were incubated with increasing concentrations of 
WhiA (indicated in /jlM above the lanes) and subjected to DNase I footprinting analysis as described in Materials and Methods. The footprint is flanked by a 
Maxam and Gilbert sequence ladder (AG). The black bracket indicates the position of the WhiA-protected region. (C) Sequence summary of DNase I 
footprinting results presented in panel B. The bracket indicates the protected sequence, and the consensus motif in shown in bold. The number indicates the 
distance (bp) between the consensus motif and the putative start codon of the whiA gene. 



tures in vitro (25). These data suggest that FilP plays a true cyto- 
skeletal role supporting tip growth. During differentiation, aerial 
hyphae arrest tip extension and initiate a massive, synchronous 
septation event, and during this switch the polarisome disappears 
(28-30). The orderly arrest of aerial growth requires whiA (5, 6), 
but the mechanisms underlying this dependence are unknown. 
One of the strongest WhiA targets, as judged by ChlP-seq, is filP 
(sven5047) (Fig. 7 and Table 1). Expression offilP is increased in a 
whiA mutant (Fig. 8), suggesting that WhiA represses filP tran- 
scription. Repression of filP by WhiA during development may 
therefore contribute to the orderly cessation of aerial growth prior 
to sporulation and thus help to explain why whiA mutants form 
abnormally long aerial hyphae. In the context of coiled-coil pro- 
teins like FilP, it is interesting that we also identified slzA (small 
leucine zipper protein A; sven5271) as a member of the WhiA 
regulon (Table 1 ). slzA encodes a coiled-coil protein of the leucine 
zipper class that interacts strongly with itself in bacterial two- 



hybrid assays, and slzA mutants show a mild developmental phe- 
notype (A. Kotun and J. McCormick, personal communication). 

Another WhiA target identified here is cslA (sven5061) (Ta- 
ble 1 ) , which encodes a putative cellulose synthase that localizes to 
hyphal tips (31, 32). CslA is proposed to reinforce growing tips 
through cellulose production, and disruption of cslA leads to a 
severe delay in the formation of aerial hyphae (31, 32). Although 
the expression of cslA is largely unchanged in a whiA mutant, there 
is a 2-fold increase at a late stage in development (20 h), suggesting 
that WhiA represses cslA (see Table SI in the supplemental mate- 
rial). cslA appears to be the first gene of a bicistronic operon that 
also includes glxA, encoding a galactose oxidase required for 
proper aerial development during osmotic stress. Although the 
two genes are cotranscribed, glxA is also expressed from an inde- 
pendent promoter within the clsA coding sequence (33). 

WhiA targets involved in cell division and chromosome seg- 
regation. Strong signals for WhiA binding were found upstream 
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FIG 7 ChlP-seq data for 12 selected WhiA target genes: ftsZ, ftsK, sepF3,ftsW,filP, whiG, cdgB, crp, sven3676, glpF, sven3229, and sven4756. Color coding of the 
ChIP samples is as follows: 3XFLAG-[Gly 4 Ser] 3 -WhiA strain during vegetative growth (V, red), at the onset of sporulation/fragmentation (F, purple), and at 
midsporulation (S, blue) and the S. venezuelae negative-control strain (C, green). Plots span approximately 3 kb of DNA sequence. Genes running left to right are 
shown in green, and genes running right to left are shown in red. The black arrow indicates the gene subject to WhiA regulation. The arrangement of the 12 panels 
mirrors that in Fig. 8. 
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FIG 8 Microarray transcriptional profiling data for 12 selected WhiA target genes {ftsZ,ftsK, sepF3, ftsW, filP, whiG, cdgB, crp, sven3676, glpF, sven3229, and 
sven4756) during submerged sporulation in wild-type S. venezuelae (purple circles) and in the congenic whiA mutant, SV11 (red triangles). In each panel, the 
x axis indicates the age of the culture in hours, and the y axis indicates the per-gene normalized transcript abundance (log 2 ). For the wild type, 10 to 14 h 
corresponds to vegetative growth, 14 to 16 h corresponds to the onset of sporulation (fragmentation), and 16 h onwards corresponds to sporulation. The 
arrangement of the 12 panels mirrors that in Fig. 7. 
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of several key genes involved in cell division in Streptomyces, in- 
cluding ftsZyftsW, and ftsK (Fig. 7 and Table 1 ). Cell division is not 
essential for viability in Streptomyces, but it is required for sporu- 
lation, which involves the synchronous formation of a ladder of 50 
or more cytokinetic rings (2). ftsZ transcription is under develop- 
mental control, being substantially upregulated from the ftsZp2 
promoter during sporulation (34). Further, in S. coelicolor at least, 
constitutive overexpression of ftsZ restores sporulation septation 
to whi mutants, suggesting that appropriate regulation of FtsZ 
levels is a vital function of the whi gene products (35). BldD was 
the first transcription factor to be identified as a direct regulator of 
ftsZ transcription, functioning to repress ftsZ expression during 
vegetative growth (14), when FtsZ is required only for the forma- 
tion of occasional vegetative crosswalk, which divide hyphae but 
do not lead to constriction and cell-cell separation. Microarray 
transcriptional profiling conducted here shows that ftsZ expres- 
sion is significantly downregulated in a whiA mutant, consistent 
with WhiA acting as a direct activator oiftsZ expression (Fig. 8). 

Similarly, our data suggest that WhiA functions as an activator 
of/fsWand/feK" transcription via direct binding to their respective 
promoters (Fig. 7 and 8; Table 1). (Note that/fsWlies in the mid- 
dle of the dew operon, suggesting that the WhiA-binding site iden- 
tifies an /fsW-specific promoter.) In Escherichia coli, FtsW is re- 
quired to recruit its cognate transpeptidase, FtsI (PBP3), to the 
division site (36), and FtsW is also likely to function as the lipid II 
"flippase" (37), transporting bacterial cell wall precursors across 
the membrane to FtsI and other transpeptidases for incorporation 
into newly synthesized septal peptidoglycan. S. coelicolor ftsW null 
mutants grown under most conditions have a white phenotype 
and make aerial hyphae that are unable to assemble FtsZ rings, 
form sporulation septa, or segregate their chromosomes (38). 
However, cell division occurs independently of/fsWand/fs/ when 
S. coelicolor is grown on low-osmolarity media (39). Since the 
S. coelicolor genome encodes four SEDS (shape, elongation, divi- 
sion, and sporulation)/penicillin binding protein (PBP) pairs, it is 
likely that another pair is able to support cell division under cer- 
tain growth conditions. Indeed, FtsZ is the only known compo- 
nent of the division machinery absolutely required for cell division 
under all conditions in Streptomyces, although it is not essential for 
viability (2, 40). 

Sporulation septa constrict over unsegregated nucleoids, and 
FtsK functions as a DNA translocase, removing trapped DNA 
from the closing septa (1-3, 41). FtsK localizes to sporulation 
septa, but it is not required for Z-ring assembly or sporulation 
septation. However, ftsK mutants have irregular DNA content in 
their spores and show greatly increased genetic instability, most 
often associated with large deletions at the ends of the linear chro- 
mosome, observations consistent with FtsK's role as a septal DNA 
pump (41). 

Our ChlP-seq data also reveal binding of WhiA upstream of 
the ftsE gene (Table 1). In E. coli, FtsE is required for constriction; 
it resembles the ATP-binding cassette of an ABC transporter, and 
it interacts with its membrane component, FtsX, as well as with 
FtsZ (42-45) . Despite the presence of a WhiA ChIP peak upstream 
of the gene,ftsE does not show strong dependence on whiA for its 
expression (see Table SI in the supplemental material). 

Chromosome segregation during sporulation also requires the 
ParAB system. ParA is a Walker A ATPase that assists ParB in the 
formation of foci formed through the binding of ParB to -20 parS 
sites that surround oriC (4, 46, 47). The observation that ParA 



forms filaments that originate from the tip and extend down the 
sporogenic cell suggests that it provides the force to organize the 
ParB-bound chromosomes prior to septation (4, 48). More re- 
cently, this process has been shown to be coordinated by a direct 
interaction between ParA and Scy, one of the components of the 
polarisome that directs tip extension in vegetative and aerial hy- 
phae (30). A previous study showed that the sporulation-specific 
upregulation of parAB transcription in S. coelicolor depends, di- 
rectly or indirectly, on several of the whi genes, including whiA 
(22), and S. coelicolor WhiA was shown to bind to the parAB pro- 
moter in vitro (9). Our microarray transcriptional profiling shows 
that the sporulation-specific upregulation of parAB transcription 
also depends on whiA and several other whi genes in S. venezuelae, 
but we did not detect binding of WhiA to the parAB promoter in 
our ChlP-seq experiments at any of the three time points exam- 
ined. svenl406 is annotated as parA-\ike, although its function has 
yet to be investigated. Our experiments identify svenl406 as a di- 
rect WhiA target and show that its expression depends on whiA 
(Table 1; see also Table SI in the supplemental material). 

SepF functions as a component of the cell division machinery 
in Bacillus, where it is suggested to play an early role in the assem- 
bly of the FtsZ ring (49), in addition to a later role in constriction 
(50). In S. venezuelae, there are three SepF-like proteins (SepFl, 
SepF2, and SepF3), but the individual (and presumably related) 
functions of these homologues are unknown. Our data show that 
the sepF3 (sven5776) promoter is a direct target of WhiA (Fig. 7 
and Table 1) and that expression of sepF3 is significantly down- 
regulated in a whiA mutant, suggesting that WhiA activates tran- 
scription of sepF3 (Fig. 8). 

Overall, these results suggest that WhiA plays a central role in 
the initiation of sporulation in Streptomyces through the activa- 
tion of multiple components of the cell division machinery. 

Other WhiA targets involved in morphological differentia- 
tion. cdgB is a target that appears to be moderately repressed by 
WhiA during the later stages of differentiation (Fig. 7 and 8; Ta- 
ble 1). cdgB encodes a diguanylate cyclase that synthesizes c-di- 
GMP, and both overexpression and deletion of cdgB inhibit the 
formation of aerial hyphae (51), the first experiments to establish 
a role for c-di-GMP signaling in Streptomyces development. Strik- 
ingly, in S. coelicolor, BldD directly regulates cdgB, in addition to 
two other genes encoding predicted diguanylate cyclases, cdgA 
and sco5511 (14, 51). The fact that cdgB is controlled by two key 
developmental regulators, WhiA and BldD, reinforces the emerg- 
ing notion that c-di-GMP signaling plays an important role in the 
regulatory network controlling differentiation in Streptomyces 
(14, 51, 52). 

ChlP-seq identified whiG (sven5300), encoding the 
sporulation-specific sigma factor o^' 6 , as a WhiA target (Fig. 7 
and Table 1). Although the full cr whlG regulon has not been de- 
fined, cr whlG has been shown to directly activate transcription of 
whiH and whil, encoding two further key sporulation transcrip- 
tion factors (Fig. 9) (53, 54). Our transcriptional profiling shows 
decreased whiG expression in an S. venezuelae whiA mutant, sug- 
gesting that WhiA activates whiG transcription (Fig. 8). Previous 
ChlP-seq studies showed that whiG is a target for regulation by 
BldD (14). Therefore, both the release of BldD-mediated repres- 
sion of whiG and WhiA-mediated activation appear to contribute 
to the regulation of whiG expression during the transition from 
aerial growth to sporulation (Fig. 9). 

It has been suggested, based on in vitro transcription, that 
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p-whiG directs transcription of the p2 promoter just upstream of 
whiA (9). However, Ainsa et al. (5) failed to detect Ea- whiG - 
directed in vitro transcription from the whiAp2 promoter under 
conditions in which the same preparation of o- whiG holoenzyme 
actively transcribed the well-validated o- wlliG target promoter, 
whilp. Further, Ainsa et al. (5) found, using SI nuclease protection 
assays, that the whiAp2 promoter is active in an S. coelicolor whiG 
deletion mutant, thus ruling out the in vivo dependence of whiAp2 
on cr whiG . It has also been suggested that WhiA and cr whiG directly 
interact and that WhiA might function as a cr whiG -specific anti- 
sigma factor (9). We have been unable to detect interaction be- 
tween S. venezuelae WhiA and cr whiG in an E. coli bacterial two- 
hybrid system. Resolution of these issues will require further 
experimental work. 

In Streptomyces, the nrdR gene encodes a repressor of the 
nrdAB and nrd] genes, which encode class la and class II ribonu- 
cleotide reductases, respectively (55-57). Ribonucleotide reducta- 
ses catalyze the reduction of nucleoside diphosphates (NDPs) into 
deoxynucleoside diphosphates (dNDPs), which are then con- 
verted into deoxynucleoside triphosphates (dNTPs), the immedi- 
ate biosynthetic precursors of DNA, and ribonucleotide reducta- 
ses play a critical role in regulating the total rate of DNA synthesis 
in vivo (58). Visualization of active replisomes using enhanced 
green fluorescent protein (EGFP) fused to a subunit of DNA poly- 
merase III (DnaN-EGFP) indicates that, during the development 
of aerial hyphae, DNA replication is dramatically upregulated in 
the long, unseptated sporogenic cell, generating the 50 to 100 
chromosomes that will eventually segregate into spores (59). This 
substantial enhancement of DNA replication presumably necessi- 
tates an increased supply of dNTPs. Here we show that WhiA 
directly binds to a region upstream of the nrdR gene and that nrdR 
expression is significantly upregulated in a whiA mutant (Table 1; 



see also Table SI in the supplemental material), suggesting that 
WhiA represses the transcription of nrdR. Decreased levels of 
NrdR would, in turn, lead to derepression of nrdAB and nrd] tran- 
scription and increased levels of the ribonucleotide reductases, 
thereby enhancing the levels of the dNTPs required for DNA syn- 
thesis in the sporogenic cell. Consistent with this logic, our mi- 
croarray analysis also reveals a corresponding indirect depen- 
dence of nrdAB and nrd] expression on WhiA prior to sporulation. 

The largest ChlP-seq peak identified in our analysis lies up- 
stream of the pyrG gene, encoding CTP synthase (Table 1), an 
enzyme that generates CTP from UTP, ATP, and glutamine. The 
expression of this gene is downregulated in a whiA mutant, sug- 
gesting that WhiA activates pyrG ( see Table S 1 in the supplemental 
material). Furthermore, we also identified pyrR (svenl088) as a 
WhiA target in ChlP-seq, and transcriptional profiling suggests 
that WhiA represses transcription of pyrR (Table 1; see also Table 
SI). In other bacteria, PyrR blocks pyrimidine biosynthesis by 
binding to mRNA derived from the pyr operon in the presence of 
uridine nucleotides (60). pyrG is not part of the pyr operon and 
therefore not regulated by PyrR. Activation of pyrG expression by 
WhiA might lead to an increase in CTP levels, but this could feed 
DNA replication in sporogenic hyphae only if it were converted to 
dCTP via three steps (CTP to CDP to dCDP to dCTP), and, apart 
from the ribonucleotide reductases, no other genes involved in 
dNTP synthesis were identified as WhiA targets. A speculative 
alternative explanation for the increase in pyrG expression during 
Streptomyces morphogenesis might be drawn from the observa- 
tion that, in several organisms, CTP synthases have the ability to 
form cytoskeletal-like filaments (reviewed in references 61 and 
62). In Caulobacter crescentus, these CTP synthase filaments ap- 
pear to have been coopted for a cytoskeletal role in determining 
cell shape, regulating the curvature of cells independently of the 
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enzyme's catalytic function (63). The mechanism by which the 
CTP synthase filaments regulate cell shape involves cross talk with 
another cytoskeletal element, the intermediate filament crescentin 
(63). It will be interesting to see whether the developmental regu- 
lation oipyrG by WhiA reflects a similar cytoskeletal role for CTP 
synthase in Streptomyces. 

Other direct WhiA targets identified here that are linked to 
morphological development include four additional genes encod- 
ing transcription factors: the cyclic AMP receptor protein (Crp), a 
PadR-family transcriptional regulator (sven3676), the SigB-like 
RNA polymerase sigma factor SigN (sven3785), and WblH 
(sven6397) (Fig. 7 and Table 1). Crp appears to function princi- 
pally as a regulator of secondary metabolism and antibiotic pro- 
duction in S. coelicolor, but several Crp targets are associated with 
morphogenesis (64). WhiA appears to repress transcription of crp 
both during vegetative growth and following differentiation 
(Fig. 8). The PadR-like regulator is likely to function as a repressor 
to influence the expression of its downstream genes in S. coelicolor, 
and both deletion and overexpression of the corresponding gene 
lead to a bald phenotype with sparse aerial hyphae (65). Tran- 
scriptional profiling in S. venezuelae reveals clear upregulation in a 
whiA mutant (Fig. 8), suggesting that WhiA represses transcrip- 
tion of the PadR-family repressor gene, cr^ 1 " 1 activity appears to be 
restricted to the "subapical stem," the cell lying immediately be- 
low the sporogenic cell in aerial hyphae, where it plays a poorly 
defined role in morphological development and the general stress 
response (66). Our ChlP-seq analysis shows that sigN is a WhiA 
target in S. venezuelae, although our transcriptomic data suggest 
that WhiA does not exert a very substantial effect on sigN expres- 
sion (see Table SI in the supplemental material). WblH is a mem- 
ber of the Wbl ( WhiB-like) family of Actinobacteria-specific pro- 
teins that also includes WhiB and WhiD (67-69). Wbl proteins 
contain four perfectly conserved cysteine residues that bind an 
oxygen-sensitive, nitric oxide (NO)-sensitive [4Fe-4S] cluster 
(70-72). The biochemical role of Wbl proteins has been a matter 
of dispute (reviewed in reference 68), but recent evidence strongly 
suggests that Wbl proteins function as transcription factors (73). 
The regulation of o- whiG , cr"^, Crp, WblH, and other transcrip- 
tion factors by WhiA emphasizes the pleiotropic potential of 
WhiA control. 

treZ, the WhiA target investigated above that lacks the GACAC 
consensus sequence but carries two copies of a GACGC motif 
(Fig. 4), encodes the metabolic enzyme TreZ, which functions in a 
pathway that breaks down the storage compound glycogen to 
yield trehalose (74). The accumulation and degradation of glyco- 
gen, stored as insoluble granules, are under developmental control 
in Streptomyces, and glycogen reserves are actively degraded dur- 
ing sporulation (75, 76), leading to an accumulation of trehalose 
that contributes substantially to the heat and desiccation resis- 
tance of mature spores (77). Thus, the activation of treZ transcrip- 
tion by WhiA (see Table SI in the supplemental material) provides 
a direct link between the regulatory cascade that controls differ- 
entiation in Streptomyces and the accumulation of trehalose in 
spores. 

WhiA targets involved in translation. Unexpectedly, several 
of the WhiA targets that we identified through ChlP-seq are in- 
volved in general translation, including the genes encoding trans- 
lation initiation factor 1 (IF-1), translation elongation factor P, 
small-subunit (SSU) ribosomal protein S8, SSU ribosomal pro- 
tein SI, and SSU ribosomal protein S10 (Table 1; see also Table SI 
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in the supplemental material). Expression of each of these targets 
is significantly downregulated in a whiA mutant, suggesting that 
an increased abundance of these components of the translation 
machinery may be a feature of morphological development in 
Streptomyces. It should be noted that SI is an RNA-binding pro- 
tein that has also been coopted as a transcription termination 
factor (78), and S10 (NusE) is involved in antitermination in 
E. coli (79). It is therefore possible that these proteins are being 
used for something other than translation in S. venezuelae. 

Conclusions. The work reported here confirms, on a genome- 
wide scale, that WhiA functions as a transcription factor and 
shows that it behaves almost equally as an activator and as a re- 
pressor. Further, supporting previous in vitro analysis (9), we 
show that most in vivo binding sites are associated with a highly 
conserved GACAC motif, but we find that there is no clear rela- 
tionship between the orientation of this asymmetric sequence and 
whether the corresponding target gene is activated or repressed by 
WhiA. 

WhiA consists of an N-terminal homing endonuclease domain 
and a C-terminal HTH domain, connected by a flexible linker (8), 
and the C-terminal HTH domain (which is completely absent 
from true homing endonucleases) binds DNA site-specifically by 
itself (9). It is therefore possible that the HTH domain of WhiA 
has entirely replaced the DNA-binding function of the homing 
endonuclease domain. However, although the homing endonu- 
clease domain cannot footprint on DNA in isolation, it does ap- 
pear to contribute to the affinity of DNA binding by full-length 
WhiA (9). Therefore, it is conceivable that the homing endonu- 
clease domain could also contribute to the specificity of DNA 
binding. It may also mediate protein-protein interaction or be a 
target for posttranslational regulation controlling the activity of 
WhiA. Extensive further research will be required to understand 
how the evolutionarily adapted structural and biochemical fea- 
tures of WhiA enable it to fulfill its bifunctional role as repressor 
and activator. 

Based on Western blot analysis, WhiA is constitutively present 
at all time points examined (Fig. 5; see also Fig. S2 in the supple- 
mental material). In contrast, based on ChlP-seq analysis, WhiA 
binds to its target promoters predominantly at the onset of spo- 
rulation during fragmentation phase (Fig. 4, 6, and 7; see also 
Fig. S3). These results suggest that WhiA activity is subject to 
posttranslational regulation. 

Counting both DNA strands, there are 15,116 GACAC se- 
quences in the S. venezuelae chromosome, but only a fraction of 
these sequences are bound by WhiA in vivo (as judged by ChlP- 
seq). In addition, there are clear targets (such as treZ) that bind 
WhiA in vivo and in vitro but do not possess this exact motif. 
Additional specificity determinants must presumably exist in ad- 
dition to the short (albeit highly conserved) GACAC motif. The 
shortness of this sequence and the consequent difficulties in bioin- 
formatic prediction of true in vivo binding sites are reminiscent of 
eukaryotic transcription factors, something that could conceiv- 
ably reflect the evolutionary history of WhiA. 

Comparing the sequences of the S. venezuelae WhiA-binding 
sites defined by DNase I footprinting with the equivalent pro- 
moter regions in S. coelicolor, the GACAC motifs within the foot- 
prints are conserved in S. coelicolor for the whiA, filP, ftsW, nrdR, 
and parA-Yikt targets but not for cslA and sven4724/4725. These 
results suggest that the WhiA regulon is likely to be broadly con- 
served between these two model species. However, because the 
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FIG 10 Comparison of the cell division machinery of E. coli/B. subtilis and Streptomyces. (A) Genes encoding key components of the divisome in E. coli and 
B. subtilis, their putative or determined function, and their presence or absence in Streptomyces. Components of the Streptomyces cell division machinery 
identified in this work as direct targets for WhiA regulation are highlighted in red and marked with an arrow. (B) Differences between the known components 
of the cell division machinery in E. coli/B. subtilis and Streptomyces. 



GACAC sequence occurs so frequently in strep tomycete genomes, 
it is not possible to make a valid global comparative analysis for 
WhiA-binding sites denned only by ChlP-seq. 

Our constructed S. venezuelae null mutant shows that WhiA is 
required for the initiation of sporulation septation and chromo- 
some segregation, as it is in S. coelicolor. It is therefore striking that 
several pivotal proteins of the Streptomyces cell division machin- 
ery, such as FtsZ, FtsW, and FtsK (Fig. 10), are encoded by genes 
that are direct targets for activation by WhiA during development. 
The initiation of the massive cell division event is tightly coordi- 
nated with the orderly arrest of apical growth in the sporogenic 
cell. It may therefore be significant that filP becomes repressed by 
WhiA at approximately the same time thatftsZ,ftsW, and/fs^Tare 
activated, given that FilP is a key component of the polarisome 
that directs polar growth at the tips of hyphae (25, 27). Further, 
recent evidence suggests that the direct interaction of another cru- 
cial component of the polarisome, Scy, with the partition protein 
ParA coordinates the arrest of tip extension in the sporogenic cell 
with the initiation of chromosome segregation prior to septation 
(30). The sporulation-specific upregulation of parAB transcrip- 
tion depends on whiA in S. venezuelae, but we did not detect bind- 
ing of WhiA to the parAB promoter in ChlP-seq, suggesting that 
this dependence is indirect. 

As illustrated in Fig. 10, many important components of the 
divisome in E. coli and Bacillus subtilis have no clear homologues 
in Streptomyces. For example, the FtsZ ring has to be anchored to 
the membrane, and in E. coli, this role is fulfilled by FtsA and ZipA 
(80), proteins absent in Streptomyces (2). This suggests that in 
Streptomyces there exist other proteins that perform functions 



central to the ability of the organism to divide and sporulate. The 
Streptomyces proteins SsgA and SsgB that function to recruit FtsZ 
during sporulation perfectly illustrate this point (81). Further cell 
division proteins unique to streptomycetes may well be under the 
control of one or more of the key developmental regulators, per- 
haps especially WhiA and WhiB, which are the only Whi regula- 
tors absolutely required for sporulation septation in S. venezuelae 
(M. J. Bibb and M. I. Buttner, unpublished data). In seeking out 
these proteins, it might be particularly fruitful to explore the over- 
lap between the regulons of WhiA and WhiB, given their sug- 
gested close interplay and the near-identical phenotypes of the 
whiA and whiB mutants in both S. coelicolor and S. venezuelae 
(Fig. 9) (5, 6; Bibb and Buttner, unpublished). 

WhiA-like proteins are found in virtually all Gram-positive 
bacteria, including many that do not sporulate, raising interesting 
questions about the evolution of WhiA and the reassortment of 
the genes under its control in different bacteria. Outside the strep- 
tomycetes, the only published genetic analysis of whiA function is 
in Bacillus subtilis. In this endospore-forming bacterium, deletion 
of whiA (yvcL) has no obvious phenotypic consequences (82), 
suggesting that WhiA does not play an important role in the con- 
trol of sporulation or cell division in this organism. It will be 
interesting to confirm that WhiA functions as a transcription fac- 
tor in Bacillus and, if so, to identify the genes it regulates. 

MATERIALS AND METHODS 

Bacterial strains, plasmids, oligonucleotides, and growth conditions. 

Strains, plasmids, and oligonucleotides used in this study are described in 
Table S2 in the supplemental material. Escherichia coli K-12 strain DH5a 
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was used for plasmid and cosmid propagation. BW25113 (83) containing 
a A RED plasmid, pIJ790, was used to create the disruption cosmid. Cos- 
mids were conjugated from the dam dent hsdS E. coli strain ET12567 
containing pUZ8002 (84) as described by Gust et al. (18, 19). Although 
ET12567 was used for this purpose, it should be noted that S. venezuelae 
does not restrict methylated DNA, and so cosmids can be conjugated from 
methylating E. coli strains. S. venezuelae strains were cultured in MYM- 
TAP (85) made with 50% tap water and supplemented with 200 ul of trace 
element solution (86) per 100 ml. Conjugations between E. coliand S. ven- 
ezuelae were essentially carried out as described for S. coelicolor ( 86 ) . How- 
ever, because S. venezuelae grows much more quickly than S. coelicolor, the 
plates were incubated either at room temperature overnight or at 30°C for 
only 7 h before overlaying with the selective antibiotics (incubation at 
30°C for longer periods before overlaying results in confluent lawns of 
S. venezuelae). Additionally, there is no need to heat shock S. venezuelae 
spores prior to the mating. Where conjugation was carried out using the 
nonsporulating \whiA::apr strain, MYM liquid cultures were set up using 
mycelium as an inoculum and actively growing cultures were concen- 
trated and used in the conjugation. 

Spores of wild-type S. venezuelae were harvested from MYM plates 
using sterile cotton pads. Three milliliters of 20% glycerol was gently 
applied across the plate using sterile forceps and the cotton pad. The 
spores were collected through the pad using a 2-ml syringe. For the non- 
sporulating whiA mutant, mycelium in 20% glycerol was homogenized 
manually 10 times in a Dounce homogenizer and appropriate (and repro- 
ducible) inocula were determined in trial experiments, as judged by opti- 
cal density (OD) and microscopic inspection. 

Construction and complementation of an S. venezuelae whiA null 
mutant. Using the Redirect PCR targeting method of Gust et al. (18, 19), 
whiA mutants were generated in which the coding region was replaced 
with a single apramycin resistance (apr) cassette. A cosmid library that 
covers >98% of the S. venezuelae genome (Bibb and Buttner, unpub- 
lished) is fully documented at http://strepdb.streptomyces.org.uk/. Cos- 
mid Sv-2-Bll was introduced into E. coli BW25113 containing pIJ790, 
and the relevant genes were replaced with the apr-oriT cassette amplified 
from pIJ773 using the primer pairs whiAdisfor and whiAdisrev (see Ta- 
ble S2 in the supplemental material). The resulting disrupted cosmids 
were confirmed by restriction digestion and by PCR analysis using the 
flanking primers whiAconfor and whiAconrev and introduced into S. ven- 
ezuelae by conjugation. Null mutant derivatives, generated by double 
crossing-over, were identified by their apramycin-resistant, kanamycin- 
sensitive, and morphological phenotypes, and their chromosomal struc- 
tures were confirmed by PCR analysis using the flanking primers whiA- 
confor and whiAconrev and by Southern hybridization using the entire 
cosmid Sv-2-Bl 1, partially digested with Sau3AI, as a probe. A represen- 
tative whiA null mutant was designated SV11. For complementation, 
whiA was amplified with the primers whiAcompfor and whiAcomprev, 
generating a fragment carrying the coding sequence and the whiA pro- 
moter, and cloned into Hindlll-cut pMS82 (87) to create pIJ6760. The 
plasmid was introduced into the whiA mutant by conjugation and fully 
complemented all aspects of the mutant phenotype. 

Construction of a 3 X FLAG-WhiA-complemented S. venezuelae 
strain. In order to engineer an S. venezuelae strain expressing a form of 
WhiA with an N-terminal, triple-FLAG tag (DYKDHDGDYKDHDI- 
DYKDDDDK), a pMS82-derived construct, pill 0600, was created via a 
two-step fusion-PCR approach. In the first step, the Sv-2-B 1 1 cosmid was 
used as a template for two separate PCRs. The first reaction amplified the 
promoter region of the whiA gene using the primer pair whiAFLAGextfor 
and whiAFLAGfusrev. The second reaction amplified the coding region of 
the whiA gene using the primer pair whiAFLAGfusfor and whiAFLAGex- 
trev. Together, the whiAFLAGfusrev and whiAFLAGfusfor primers con- 
tain the sequence encoding the triple-FLAG tag via a 24-bp overlapping 
section. In the second step, a PCR using the nested primers whiAFLAG- 
nesfor and WhiAFLAGnesrev was used to amplify the entire whiA gene 
and its promoter, fusing the two products from step 1 together and 



incorporating the 3 X FLAG tag sequence between them. The nested 
primers whiAFLAGnesfor and WhiAFLAGnesrev additionally contain 
the Hindlll and Kpnl sites, respectively, to enable cloning into Hindlll-, 
KpnI-cut pMS82. The resulting vector was named pI110600. In order to 
insert a [Gly 4 Ser] 3 linker between the 3 X FLAG tag and the coding region 
of whiA, a second fusion experiment was conducted using pIJ10600 as a 
template. Once again, two PCRs were carried out in a first step using the 
primer pairs whiAFLAGnesfor/whiALINKfusrev and WhiAFLAGnesrev/ 
whiALINKfusfor. Together, the whiALINKfusrev and whiALINKfusfor 
primers contain the sequence encoding the [Gly 4 Ser] 3 linker (GGTGGA 
GGCGGTTCAGGCGGAGGTGGCTCTGGCGGTGGCGGTAGT) via a 
24-bp overlapping section. The second step, employing the same nested 
primers, resulted in fusion of the two products in the first step and incor- 
porated the linker sequence between them. Restriction digestion using 
Hindlll and Kpnl, followed by ligation into Hindlll-, KpnI-cut pMS82, 
created pIJ10601. The plasmids pIJ10600 and pIJ10601 were introduced 
individually into the Aw/nA mutant by conjugation, and their ability to 
restore sporulation was assessed both on solid and in liquid MYM. 

RNA isolation and DNA microarray analysis. RNA isolation from 
S. venezuelae and DNA microarray analysis were performed as described 
previously (16, 88). The CEL files received from the scanner were read into 
the R package for statistical computing (89) using the ReadAffy function 
of the affy package (90). The rma function of the affy package was used to 
compute an ExpressionSet object from the CEL files. This ExpressionSet 
object contains the expression values (log 2 ) for each gene in each CEL file. 
The function ImEit of the limma package ( 9 1 ) , along with a suitable design 
matrix, was used to combine replicate arrays into single coefficients of 
expressions for each gene at each time point or strain into an MArrayLM 
object. Expression values were retrieved from the MArrayLM object and 
subjected to a per-gene normalization to the median before being used to 
generate the graphs shown in this paper. 

Chromatin immunoprecipitation. S. venezuelae ATCC 10712 and 
S. venezuelae SVll-pIJ10601 (AwhiA::apr flttB Mrj ::3XFLAG- 
[Gly 4 Ser] 3 -vf/nA) were grown in four 30-ml volumes of MYM-TAP for 
the appropriate length of time. Formaldehyde was added to cultures at a 
final concentration of 1% (vol/vol), and incubation was continued for 
30 min. Glycine was then added to a final concentration of 125 mM to stop 
the cross-linking. The samples were left at room temperature for 5 min 
and washed twice in phosphate-buffered saline (PBS) buffer (pH 7.4). The 
pellets were resuspended in 0.5 ml of lysis buffer (10 mM Tris-Cl, pH 8.0, 
50 mM NaCl), containing 15 mg ml -1 of lysozyme and IX protease 
inhibitor (Roche Applied Science) and incubated at 37°C for 25 min or 
until lysed. Subsequently, 0.5 ml immunoprecipitation (IP) buffer 
(50 mM Tris-Cl, pH 8.0, 250 mM NaCl, 0.8% [vol/vol] Triton X-100), 
containing IX protease inhibitor, was added, and the samples were 
chilled on ice. The samples were sonicated for six cycles of 20 s each at 
10 u,m to shear chromosomal DNA into fragments ranging from 300 to 
1,000 bp on average. The samples were centrifuged twice at 16,000 X g at 
4°C for 10 min to clear the cell extracts, after which 50 ul of each was set 
aside for total-DNA extraction. The remaining lysates were added to M2 
(Sigma-Aldrich A2220) gel suspension (40 ul per ml of lysate), prepared 
and washed as described in the manufacturer's instructions. The mixtures 
were then incubated on a rotating wheel at 4°C overnight. The samples 
were centrifuged for 30 s at 4°C and 5,500 X g, and the pellets were washed 
twice with 0.5 X IP buffer and then twice with IX IP buffer and trans- 
ferred to new Eppendorf tubes after the first washing step. The pellets and 
50 ul of total cell extracts (set aside earlier) were eluted overnight at 65°C 
in 100 ul of IP elution buffer (50 mM Tris-Cl, pH 7.6, 10 mM EDTA, 1% 
SDS) to reverse the cross-links. The samples were centrifuged at 16,000 X 
g for 5 min to remove the beads. The pellets were reextracted with 50 ul of 
TE buffer (10 mM Tris-Cl, pH 7.4, 1 mM EDTA) and incubated with 
0.2 mg ml -1 proteinase K (Roche) for 2 h at 55°C. The samples were 
extracted twice with phenol-chloroform and once with chloroform and 
further purified using QiaQuick columns (Qiagen). DNA was eluted in 
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50 ul EB buffer (10 mM Tris-Cl, pH 8.5) and quantified using a Nano- 
Drop spectrophotometer (Thermo Scientific). 

Library construction and sequencing. Library construction and se- 
quencing were performed by The Genome Analysis Centre (TGAC), Nor- 
wich Research Park, Norwich, United Kingdom. The TruSeq ChIP sample 
preparation kit from Illumina Inc. (part 15023092, revision A, August 
2012) was used according to the recommended protocol. In brief, 5 to 
10 ng ChIP DNA is blunt ended and phosphorylated, and a single "A" 
nucleotide is added to the 3' ends of the fragments in preparation for 
ligation to an adapter with a single-base "T" overhang. The ligation prod- 
ucts are purified and accurately size selected by agarose gel electrophore- 
sis. Size-selected DNA is then purified and PCR amplified to enrich for 
fragments that have adapters on both ends. The final purified product is 
then quantified using a combination of Bioanalyzer DNA HS Chip from 
Agilent Technology 2100 Bioanalyzer and Qubit 2.0 from Invitrogen. The 
libraries were normalized to 10 nM and pooled, and quantitative PCR 
(q-PCR) was performed prior to cluster generation. 

Each library pool was diluted to 2 nM with NaOH, and 5 fjl was 
transferred into 995 fjl HT-1 to give a final concentration of 10 pM. The 
normalized library (120 /oil) was then transferred into a 200-/J.1 strip tube 
and placed on ice before loading onto the Illumina cBot. Flow cells were 
clustered using TruSeq Single-Read Cluster Generation Kit v3, using the 
SR_Amp_Lin_Block_Hyb_v8.0 recipe. Following the clustering proce- 
dure, the flow cell was loaded onto the Illumina HiSeq2000 instrument 
according to the manufacturer's instructions. The sequencing chemistry 
used was TruSeq SBS version 3 using Illumina software HCS 1.5.15.1 and 
RTA 1/13/48. Each library pool was run on the Illumina HiSeq2000 in- 
strument with 100-bp single-read run metrics. 

ChlP-seq data analysis. The reads in the fastq files received from the 
sequencing contractor were aligned to the S. venezuelae genome (Gen- 
Bank accession number NC_018750) using the bowtie2 software (92), 
which resulted in one SAM (.sam) file for each fastq file. All further oper- 
ations described below were carried out using a combination of Perl 
scripts dependent on the BioPerl toolkit (93) and R (89) scripts. From 
each SAM file, coverage at (number of reads mapping to) each nucleotide 
position of the S. venezuelae genome was calculated and the output was 
saved in files here referred to as coverage files. For each coverage file, a 
local enrichment was calculated in a moving window of 51 nucleotides 
(nt) moving in steps of 25 nucleotides as (the sum of coverage at each 
nucleotide position in the 5 1 -nt window) divided by (the sum of coverage 
at each nucleotide position in a 2,001 -nucleotide window centered 
around the 51 -nucleotide window). This results in an enrichment ratio 
value for every 25 nucleotides along the genome. The enrichment ratios 
thus calculated were stored in files referred to as enrichment files. 

In each difference file, the P value of each difference was calculated 
assuming that the differences were normally distributed using the pnorm 
function of R. 

Now, for each pair of total DNA and IP DNA enrichment ratios, the 
values of total DNA were subtracted from the values of IP DNA and the 
differences stored in files here referred to as the difference files. Thus, we 
arrived at a difference file for each of the three time points for the FLAG- 
tagged pairs and one file for the wild-type (WT) pair. In each difference 
file, the P value of each difference was calculated assuming that the differ- 
ences were normally distributed using the pnorm function of R. Differ- 
ences with adjusted P values (using the pnorm function of R with the 
Hochberg method) less than or equal to le— 4 were saved to new files 
referred to as sigDiff files. All the four types of files, i.e., coverage, enrich- 
ment, difference, and sigDiff, are two-column files where the first column 
is always position on the genome and the second column contains the 
relevant value. 

The gene coordinates in the GenBank file and the genomic positions in 
the sigDiff files were combined to arrive at lists of genes on the left and 
right of the significantly different genomic positions. The relationships of 
the difference values and these genes were manually verified by viewing 
the difference files and the features in the S. venezuelae GenBank file as 



tracks in the Integrated Genome Browser (94). Significant different 
genomic positions located more than 300 bp upstream of the nearest start 
codon were separated from the data set. The most significant genomic 
position for each transcription unit was then selected. Finally, the infor- 
mation above was combined with Affymetrix microarray data to generate 
Table S 1 in the supplemental material. 

Western blotting, DNase I footprinting, and microscopy. For details 
on Western blotting, DNase I footprinting, and microscopy, see Text S 1 in 
the supplemental material. 

SUPPLEMENTAL MATERIAL 
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